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Abstract 

The production of electron-positron pairs in a vacuum neutron star magnetosphere is investigated 
for both low (compared to the Schwinger one) and high magnetic fields. The case of a strong 
longitudinal electric field where the produced electrons and positrons acquire a stationary Lorentz 
factor in a short time is considered. The source of electron-positron pairs has been calculated 
with allowance made for the pair production by curvature and synchrotron photons. Synchrotron 
photons are shown to make a major contribution to the total pair production rate in a weak 
magnetic field. At the same time, the contribution from bremsstrahlung photons may be neglected. 
The existence of a time delay due to the finiteness of the electron and positron acceleration time 
leads to a great reduction in the electron-positron plasma generation rate compared to the case 
of a zero time delay. The effective local source of electron-positron pairs has been constructed. It 
can be used in the hydrodynamic equations that describe the development of a cascade after the 
absorption of a photon from the cosmic gamma-ray background in a neutron star magnetosphere. 
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1. INTRODUCTION 



Over more than 40-year-long history of studying radio pulsars, attention has been focused 
on the investigation of normal pulsars. However, there exist nonstationary radio sources 
associated with neutron stars. The characteristics of their radio emission vary with time, 
while the radio emission itself can disappear for a while. In particular, intermittent pulsars 
and rotating radio transients can be classified as nonstationary radio sources. Observational 
data on the first intermittent pulsar were pubhshed in 2006 [1]. A significant detected effect 
is the difference between the neutron star spin-down rates in the periods of emission and 
silence. Analysis of archival data from the Multibeam Pulsar Survey with the Austrahan 
Parkes radio telescope has revealed rotating radio transients (RRATs) — the sources of short 
and relatively bright single radio bursts [2]. The width of individual bursts is 2 — 30 ms and 
the fiux density reaches 10 Jy at 1.4 GHz [3]. For most radio pulsars, the fiux density at 
this frequency does not exceed several mJy (1 Jy = 10~^^ Wm~^ Hz~^). Note that RRATs 
can also be observed at a low radio frequency of 111 MHz [4]. Although the time interval 
between two successive bursts is random, it is a multiple of some fixed period that lies within 
the range 0.1 — 6.7 s and is believed to be the neutron star rotation period. This rotation 
period exceeds that of typical radio pulsars, which is ~ 0.3 s. 

Generally, it is rather difficult to detect nonstationary radio sources and to determine 
their characteristics. One of the main reasons for this is their long stay in a state of silence. 
Such silence is also observed in nulling pulsars; the ratio of the nulling time to the total time 
of pulsar observation can reach 95% [5]. However, silence itself is not a necessary condition 
for nonstationarity, because radio emission nonstationarity also manifests itself in ordinary 
radio pulsars as mode switching. Observations in recent years suggest that the spin-down 
rate of a pulsar also changes in the absence of its switch-off [6]. This change is abrupt, is 
often quasi-periodic, and correlates with the observed change in pulse shape for six pulsars. 

Magnetars arc also characterized by a transient radio emission. Magnetars arc neutron 
stars with superstrong surface magnetic fields of the order of or higher than the Schwinger 
magnetic field. Their activity manifests itself mainly in bright X-ray and gamma-ray bursts. 
Their radio emission is observed at both low [7-10] and high [11, 12] radio frequencies. The 
radio emission could be related to X-ray bursts, but there exists a magnetar, PSR J1622- 
4950, without any observable bursty behavior in the X-ray band [13]. Moreover, based on 
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archival data, we can assert that this radio source can be switched off for a time of about 
several hundred days. Since plasma outflows from the magnetosphere are currently believed 
to be responsible for the observed radio emission from neutron stars (see, e.g., the review 
[14]), the activity of the mentioned radio sources suggests a possible cessation of the plasma 
generation in their magnetospheres [15]. Analysis of observations for the intermittent radio 
pulsars PSR B1931+24 and PSR J1832+0029 shows that the observed change in neutron 
star spin-down rate in the absence of radio emission can be explained only if the neutron 
star magnetosphere becomes a vacuum one: the plasma not only is not generated in the 
magnetosphere but also does not come from the neutron star surface [15]. However, the 
plasma production through the absorption of photons from the external cosmic gamma-ray 
background can occur in a vacuum magnetosphere [16-18]. 

Let there be a primary galactic photon whose energy and propagation direction are such 
that the transverse momentum component exceeds 2meC, where mg is the electron mass and 
c is the speed of light. This photon can then produce an electron-positron pair in a magnetic 
field. In general, the pair particles are produced at high Landau levels and, passing to the 
zeroth Landau level, emit synchrotron photons. In a vacuum magnetosphere, there is a 
strong longitudinal electric field E^y whose characteristic strength is defined by the relation 

E\\/B ~ VlRs/c = Rs/Rl ~ 10"^ 

where Rs ~ 10 km is the neutron star radius, Rl = c/Q is the light cylinder radius, and Q 
is the angular frequency of neutron star rotation. The particles in such a field acquire some 
stationary Lorentz factor 79. This Lorentz factor can be determined from the condition of 
balance between the work of the longitudinal electric field and the curvature radiation power 
and can reach 70 ~ 10^ — 10^ in a vacuum field [16, 17]. We have curvature and synchrotron 
photons that, in turn, produce electron-positron pairs. The secondary particles also emit 
synchrotron photons, acquire the stationary Lorentz factor 70, and emit curvature photons, 
which the particles of the primary pair continue to do. Subsequently, the next generation of 
particles is produced and so on. To study the development of an electron-positron plasma 
generation cascade triggered by the absorption of a photon from the external cosmic gamma- 
ray background, it is necessary to calculate the source of electron-positron pairs. This is 
the goal of our work. Here, it should be pointed out that the problem under consideration 
differs significantly from the problem of stationary electron-positron plasma generation in 
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a neutron star magnetosphere filled with a dense plasma [19, 20]. A strong longitudinal 
electric field in the latter case exists only in the polar region of the magnetosphere near the 
neutron star surface, so that the plasma is produced only in this bounded region and the 

generation cascade is usually limited to two generations. 

The paper is structured as follows. In Section 2, we derive the integral equation for a self- 
consistent source of electron-positron pairs by taking into account the particle production by 
curvature and synchrotron photons and find its solution. In Section 3, we consider the initial 
stage of development of the pair production cascade when the hydrodynamic equations are 
not yet applicable. In Section 4, we investigate the influence of the time delay between the 
production of particles and the absorption of the photons emitted by them in a magnetic 
field on the electron-positron plasma generation rate and calculate the effective local source 
of electron-positron pairs. Our main conclusions are presented in Conclusions. 



2. PAIR PRODUCTION 

Consider the filling of the magnetosphere of a neutron star with a surface magnetic field 
B ~ 10^^ G. We assume the initial state of the magnetosphere to be a vacuum one. This 
implies that the existing charge density in the magnetosphere is small compared to the 
Goldreich- Julian density poj — ■ B/27rc, when the longitudinal electric field is screened. 
The corresponding electromagnetic field was calculated in [21]. If we assume that there is 
no free charge escape from the stellar surface, then the magnetosphere can be filled only 
through the production of electron-positron pairs by high-energy gamma-ray photons from 
the external Galactic background. Allowance for the influence of these photons on the 
particle production processes in pulsars was shown to be necessary in [22]. It was pointed 
out in the latter paper that an influx of new charges that can be provided by the absorption 
of photons from the external cosmic gamma-ray background is needed in the Ruderman- 
Sutherland model [23] to trigger each new spark in the polar gap of a pulsar magnetosphere. 
We will add that when the polar cap of a radio pulsar is considered, the case of free particle 
escape from the surface of a neutron star with a nondipolar magnetic fleld is also investigated 
[24, 25]. 

In what follows, we will measure the magnetic fleld strength in units of the Schwinger 
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(critical) field 

9 

Bcr = ~ 4.414 X 10^^ G, 

en 

where e is the positron charge. As the units of mass, length, and time, we will take, 
respectively, the electron mass iJie ~ 9.109 x 10~^® g, the Compton electron wavelength 
X = h/nieC fa 3.862 x 10"^^ cm, and its ratio to the speed of light X/c 1.288 x 10~^^ s. 
Formally, this means that H — X — c — 1. Thus, the length, velocity, time, and energy will 
be measured in units of X, c, X/c, and mgC^, respectively. We will note at once that in these 
units 1 cm ^ 2.590 x 10^° and 1 s ^ 7.763 x lO^o. 

Wc consider the production of an electron-positron plasma in a weak magnetic field 
5 -C 1. In this case, the number of electron-positron pairs produced per unit time per unit 
volume in a unit interval of longitudinal Lorentz factors is given by the expression [26] 

2a / a\ , . 

where qphik) is the number of photons produced per unit time per unit volume in a unit 
energy interval, a = 4/3B ~ 100, 7|| is the longitudinal Lorentz factor of a particle, and 
A ~ 10 ^ 1 is the logarithmic factor. We see that the source of particles Qp is proportional 
to the source of photons Qph, which is represented as the sum of the source of curvature 
photons Qcurv s-nd the source of synchrotron photons Qsyn' 

Qph{k) = Qcurvik) + Qsynik). (2) 

The curvature photons are produced by ultrarelativistic particles as they move along 
curved magnetic field lines and have a characteristic energy 

k - (3) 
z p 

where p is the radius of curvature of the particle trajectory, which virtually coincides with the 
radius of curvature of the magnetic field lines. The number of curvature photons produced 
per unit time in a unit energy interval near k is given by the expression 

Pcurv{lf\\,k) = ^ ( ) ' 

V o TT 7|| \ f^curv / 

where 

oo 

(p{x)= / K5/3{y) dy, (5) 



-^5/3 (y) is the 5/3-order Macdonald function, a. — e^/hc ^ 1/137 is the fine-structure 
constant. 

Let us introduce the spectral particle distribution function F{j\\) equal to the number 
of particles per unit volume per unit interval of longitudinal Lorentz factors. The source of 
curvature photons can then be represented as 

oo 

Qcurvik)^ J Pcurv{lhk)F{-f\\)d-f\l, (6) 

Tmin 

where 7inin ~ 100 is the minimum Lorentz factor of the produced particles [26]. 

Let us calculate the source of synchrotron photons qsyn{k)- Let initially there be a 
photon with an energy ki propagating at an angle x to the magnetic field direction. An 
electron-positron pair whose particles initially have a Lorentz factor 7, and move at an 
angle Oi to the magnetic field direction is produced after its absorption. Subsequently, the 
particles lose their transverse velocity component through synchrotron radiation and pass 
to the zeroth Landau level, while having some final Lorentz factor 7/ = 7|| that we call 
the longitudinal one. The law of conservation of the energy 7, = ki/2 [27] and the law of 
conservation of the longitudinal momentum component fcj cos x = where v\\ = Vi cos 9i 

is the longitudinal electron or positron velocity component and Vi is the initial particle 
velocity, hold at the instant of electron-positron pair production. It is important to note that 
the longitudinal velocity component remains unchanged during the emission of synchrotron 
photons: I'll = vcos9, where v — ^/l — 1/7^ and 9 are some intermediate values of the 
particle velocity and pitch angle, 7 is an intermediate Lorentz factor of the particle. The 
photons propagate at small angles to the magnetic field, x ^ 1. Consequently, the pitch 
angles of the particles produced by them arc also small, because ^ 9 ^ 9i < x- Using 
the fact that, in this case, sin ^ w ^, we have the following dependence of the particle pitch 
angle on its current Lorentz factor: 

Setting the pitch angle equal to zero gives the longitudinal Lorentz factor 

711 = 1/X- (8) 
To simplify the calculations, it is convenient to introduce an angle rj via the relation cos rj — 
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1/7X- As we see from Eq. (7), the pitch angle is then given by the expression 

^ = Xsin77, (9) 

with r]i = aiccos{2 / kix) corresponding to the initial pitch angle 6i. 

The number of synchrotron photons emitted by a particle with a Lorentz factor 7 and a 
pitch angle 9 per unit time in a unit energy interval near k is 

P,,„(^,7,^) = ^4<^fAV (10) 



v/37r7^ \ksyn/ 

where 

3 

ksyn = 2^^^^ 

is the characteristic energy of synchrotron photons. During the emission of synchrotron 
photons, the particle energy decreases: 

'i = -la(BO,f. (U) 

Here, the total intensity of the synchrotron radiation obtained by integrating Eq. (10) with 
the weight k over all energies from to 00 appears on the right-hand side. 

To find the total number of photons Afsyn in a unit energy interval near k emitted by 
one particle in the synchrotron radiation time, Eq. (10) should be integrated over all times 
by taking into account the time dependence of the particle Lorentz factor and pitch angle. 
Integration over the angle rj from rji to corresponds to integration over the time from 
to 00. Given Eqs. (7), (9), and (11), it is easy to obtain the expression 

m 

J^syn [h: X, k) = / dri— (p a— — : . 

21: J sm 77 \ 2 sm 77 / 


The spectral distribution Msyn depends on the energy ki of the initial photon that pro- 
duced an electron-positron pair and its pitch angle x- However, these two quantities are not 
independent and their relation is given by the condition for photon absorption in a magnetic 
field [26]: 

kiX^ a ^ 1 
2 A cos rji 

Since the relation a/A ~ 10 » 1 holds, the angle rji is close to 7r/2 or, more specifically, 

TT A 

2 a 
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Given Eq. (8), the source of synchrotron photons can be written as 

oo 

Qsynik) ^2 J qp{-f\\)Afsyn (^^TUp dj\\. (12) 

Tmin 

The electron-positron pair production rate Qp is proportional to the total source of photons 
Qph containing the part responsible for the synchrotron photons. In turn, the synchrotron 
photon production rate qsyn itself is proportional to the pair production rate qp. Given 
Eqs. (1), (2), and (12), this self-consistency condition allows the integral equation to find 
the source of electron-positron pairs Qp to be derived: 

oo rii 

%(7||) = x^c...(27|,^j+^ey ^^^Ay— ^j' (13) 

Tmin 

where 

9x/3 
^ 87r ■ 

Consider some volume of an electron-positron plasma in a magnetic field. The plasma 
production process can be represented as follows. Let there be some photon, a curvature or 
synchrotron one, that produces an electron-positron pair. After the emission of synchrotron 
photons, the particles at the zeroth Landau level have some longitudinal Lorentz factor 
and, consequently, emit curvature photons. As we see from Eq. (6), the source of curvature 
photons depends on the spectral distribution of particles in longitudinal Lorentz factor. 
Below, we will consider the case where a fairly strong longitudinal electric field exists. The 
electron and the positron then accelerate in the time 

7o 

and move in opposite directions with the characteristic Lorentz factors 

(see [16, 17]). Here, we do not assume that the field is equal to the external electric field 
in a vacuum magnetosphere. This field is the total longitudinal electric field that is the sum 
of the external field and the plasma screening field. The field strength is assumed to be 
sufficient for electron-positron pairs to be efficiently produced. Thus, we can write the total 
particle distribution function as 

F(7||) = (n+ + n_)5(7|| - 70), (14) 



where n+ and n_ are the positron and electron number densities, respectively, and 5{x) is 
the delta function. This by no means implies that the source of particles qp{'j\\) must be 
proportional to the delta function, because it depends on the source of photons without 
a monochromatic spectrum. Nevertheless, the produced particles will acquire the Lorentz 
factor 7o irrespective of the initial particle spectrum. For this reason, it is important for us 
to find the total number of electron-positron pairs produced per unit time per unit volume. 

The particle distribution function (14) is proportional to the total particle number den- 
sity. Since the integral equation (13) is linear, the source of electron-positron pairs is also 
proportional to n+ -|- n_ . Below, it is convenient to assign the function to one particle by 
formally assuming in all of the previous formulas that n_|_ +n_ = 1. In particular, the source 
of curvature photons Qcurvik) will be given by Eq. (4), in which we should set 7|| — Jq. The 
number of electron-positron pairs produced per unit time per particle is then 

oo 

<9(7min) = J gp(7||)c?7||- (15) 



7i 



min 



Above, we have not mentioned that, apart from curvature and synchrotron photons, the 
produced particles also emit bremsstrahlung photons because of their acceleration in the 
longitudinal electric field . The bremsstrahlung intensity is [28] 

2 2 

It should be compared with the intensity of the curvature radiation 

2 7^ 

where 7 is the particle Lorentz factor. When the condition of quasi-stationary motion 
Wcurv — -E'li corresponding to 7 = 70 is met, we have 

= laE^^ « 1. (17) 

At £'11 ~ 10~^ this ratio is equal in order of magnitude to Whrems/Wcurv ~ 10^^. Obviously, 
condition (17) is always met, because the longitudinal electric field is definitely lower than 
the Schwinger one, E\\ < 1. Otherwise, we would have the direct production of electron- 
positron pairs from a vacuum, which would immediately lead to electric field screening. 

However, the particle initially has a Lorentz factor much smaller than 70. The brems- 
strahlung intensity will then still exceed the intensity of the curvature radiation until the 



Lorentz factor reaches ■\/pE\\. For E\\ ~ 10"^ and p ~ 10^^, the characteristic transition 
Lorentz factor is 10^. To estimate the total contribution from bremsstrahlung photons, let 
us find the total bremsstrahlung energy over the entire time of particle acceleration to the 
stationary Lorentz factor 70. In this case, we assume that the initial particle Lorentz factor is 
much smaller than 70. This will allow us to estimate an upper limit for the bremsstrahlung 
energy, because, obviously, it is proportional to the difference between the final and ini- 
tial Lorentz factors of the particle being accelerated. The total energy of the emitted 
bremsstrahlung photons in the time Tst will be 

The total energy of the curvature radiation in the particle acceleration time is obtained by 
integrating Eq. (16), in which the current Lorentz factor 7|| — E^^t should be taken instead 
of 7, over the time from to Tgt' 

■^curv ''curv - ■ 
O 

Now, it is easy to obtain the relation 

Eljrems r ^^brems t 

= 5 <^ 1. 

Ecurv ^^curv 

Thus, the ratio of the total energy of the bremsstrahlung photons to the total energy of 
the curvature photons emitted in the time of particle acceleration to the stationary Lorentz 
factor 7o is very small and does not depend on the value of 70 itself. We will add that there 
is an additional contribution from synchrotron photons to the total pair production rate in 

a relatively weak field B 1, which, as a rule, exceeds the contribution from curvature 
photons. Consequently, the presence of bremsstrahlung photons may be disregarded. 



2.1. The Mellin Transform 

The integral equation (13) can be solved as follows. Consider the second term on the 
right-hand side of Eq. (13). In this double integral, we will change the order of integration 
and initially consider the integral 



00 




'Ymin 
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where 

^ 0? cos^ ri 
T = -^7||— ^- 

A SUIT] 

The quantity T is hmited from below by 

T ^ A711 ^ A7„,i, > 1 

and, in any case, T ^ 1000. The integral / is a repeated one (see (5)) and, once the order 
of integration has been changed, can be represented as 

I = ) qipi^rain) + dy Kr,i^{y) gi^, I — J , (18) 

\ 7min / J \y / 



where we introduced a function 

qip{x) 



00 



Below, we will consider the functions qp{x) and qip{x) for all values of the argument x from 
to 00. We do not restrict ourselves to x ^ 7inin at which the function qp{x) is initially defined 
and assume that qp{x) was somehow extended to the interval (0, 7niin)- For us, it is only 
important that the function qp{x) coincides with the solution of the integral equation (13) 
for a; = 7|| ^ 7min- 

Let us take a power function as a trial one: 

qip{x) = x~\ (19) 

Here, s is some complex number such that integral (18) is a priori assumed to be convergent. 
Let us introduce the power moments of the Macdonald function K^j3{x) using the relations 



00 

Ms{x) = J y'K,/3{y)dy, 

X 

M,.M,(o) = 2-r(| + i)r(|-i 



3^ 

where r{x) is Euler's gamma function [29]. The integral on the right-hand side of Eq. (18) 
will then be 

' T 



J dy K,/s{y) qip (^^^ = T"^ (^M, - M, 



7n 
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The relation 

^ A> 1 

Tmin 

allows us to simplify the calculations and to use the asymptotic representation of the Mac- 
donald function at large values of the argument [30] 

/ TT _ 

-^^^5/3(2;) Y ^ ^ 00. 

Using this relation, we obtain an asymptotic expression for the incomplete moments of the 
Macdonald function and, in particular, for the function (p{x): 

TT 



Ms{x) ~ ^-x'-^/^e-"", x^oo, (20) 

ip{x) = Mo{x) ~ y^e-", x^oo. (21) 
Using Eqs. (18)-(21), we then obtain 

or, representing / as a linear integral operator I[qip{x)]{'y\\) (see Eq. (18)), 

I[x-^]{y)^L{s)y-^, (22) 

where 

/^cos^,\- 
^ ^ \A sin 77 J 

Let there be some function f{x) defined in the interval (0, 00). This function can be 
associated with a function 

00 

g{s) = mfi^m = J x^-'f{x) dx, (23) 



where s — a + it is a complex parameter. The integral transform (23) is called the Mellin 
transform of the function f{x). If the transform g(s) is known, then the inverse transform 
f{x) can be found from the formula for the inverse Mellin transform: 

(T+ioo 

f{x) = m-'[gis)]ix) = ^^ J x-'g{s)ds, (24) 
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where cr = Res and x > 0. The function f{x) must satisfy some very weak constraints. 
The following constraints can be taken as the latter [31]. Let f{x) be absolutely integrable 
in any finite interval {£i,£2), < £i < £2 < oo, and satisfy the constraints < Ax~'^^ 

for < a; ^ £1 and < Ax~'^'^ for x ^ 62, Ui < o"2, A = const > 0. The transform g{s) 

then exists and is an analytic function in the band o"i < Res < o"2. Formula (24) is valid at 
each continuity point of the function f{x). The integration in Eq. (24) is over any vertical 
straight line in the analyticity band and the integral at infinity is understood in the sense 
of a principal value. Note that there also exist generalizations of the integral transform (23) 
to the space of generalized functions [32]. 

We can now consider the case of arbitrary functions qp{x) and qip{x) by assuming that 

qp{x) = 9Jt-^[M(s)](x), qip{x) = m-\Mi{s)]{x). 

The sought-for functions A^(s) and M.i{s) are related by the obvious relation 

Mi[s) = — —. 

s 

The integral equation (13), once it has been represented as 

Vi 

M-'[M{s)]{^lO = ^q^ur. (2711^) + dv'-^ I [m-'[Mi{s)]{x)] (7||) (25) 



can be easily solved. The hnearity of the operator I and relation (22) allow the following 
expression to be derived: 

/ [m-'[Mi{s)]{x)] (7||) = 9Jt-^[L(s)Mi(s)](7||). (26) 

Now, it is necessary to take into account Eq. (26) and then apply the direct Mellin transform 
operator 9Jt to both sides of Eq. (25). Using the formula 

mMx)]{s) = ^ (27) 

S 

and the property 

m[f{Ax)]{s) = A-m[f{x)]{s), A>0, (28) 

when calculating the Mellin transform of the source of curvature photons qcurv, we ultimately 
obtain the expression for the Mellin transform A4{s) of the function qp{x) 

13 



Here, ko — ?>^l/2p is the characteristic energy of the curvature photons (see (3)) correspond- 
ing to 7|| = 7o and the function Vg is defined by the relation 

Vi 







2 \ 1— •s 

cos 77 ^ 



smT^ 



It was shown in [33] that M.{s) is a meromorphic function and, consequently, has no other 
singular points, except the poles, in the complex plane. To choose the path of integration 
in Eq. (24), it is necessary to find the analyticity band of the function M.[s). Obviously, 
this band must be bounded on the left and the right by the vertical straight lines on which 
the function has singular points, which are poles in the case under consideration. We will 
search for the bounding poles on the real axis. In addition, we will be interested in s > 2/3. 
The poles will then be specified by the equation 

P{s) = 1, (30) 

where we introduced a function 



Let Smin and Sma.x be the left and right bounding poles, respectively. The cumbersome 
calculations performed in [33] allow us to find that 

_ A 

"Smin — 2 - —, (31) 

am A 



'max 



= 2 + eA. (32) 

We see that the left bounding pole is always leftward of s = 2. The quantity A/aln A 0.05. 
For the right bounding pole, we have s^ax 1- At A ~ 10, s^ax ~ 30. The estimated upper 
hmit for the absolute value of the function P{s) leads us to conclude that \P{s)\ < 1 inside 
the band Smm < Res < s^ax [33]. Since relation (30) holds at the poles, the function 
M.(s) inside the same band has no poles. Thus, the function Ai{s) is regular in the band 
Smin < Re s < Smax SLiid the path of integration in Eq. (24) can be taken in the form of a 
vertical straight line parallel to the imaginary axis and located inside this band. 

Before we calculate the function (5(7min), let us investigate the properties of the produced- 
particle spectrum qp{x). We can consider the behavior of this function at low and high 
values of the argument. We will be interested primarily in the low-energy asymptotics of 
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the spectrum. In order to have the right to consider it, we will assume the existence of some 
transition Lorentz factor ^tr that separates the entire spectrum into two regions: the region 
of low energies and the region of high energies. Since the entire spectrum extends from 7 ^;^ 
to CX3, for such a separation to be possible, we must assume that at least 7min < 7tr- 

Let us calculate the low-energy asymptotics of qp{x). For this purpose, it is convenient 
to express the Mellin transform 9Jt in terms of the equivalent two-sided Laplace transform 
2, specified by the formula 

00 

g{s)=2[f{x)]{s)= j e~^'f{x)dx. 
—00 

The corresponding inverse two-sided Laplace transform is 

a+ioo 

fix) = 2-'[g{s)]ix) = ^ / e^^gis) ds. (33) 

a—ioo 

Then, 

qp{x)^2-^[M{s)]{-\nx). 

Consider the range of low x. lix ^ 0, then — In a; > and, moreover, —Inx^l. As follows 
from Eq. (33), it is advantageous to have s such that Re(s) are as small as possible. This can 
be achieved by transferring the path of integration as far to the left as possible. Obviously, 
the integral is invariable in this case if the path is transferred within the regularity band. 
To obtain the principal term in the asymptotic expansion of the function qp{x) at a; = 0, the 
path should be transferred through the left bounding pole Smin- The original integral will 
then be equal to the sum of the residue of the integrand at s = Smin multiplied by 27ri and 
the integral over the transferred path. The former term will give the principal term of the 
asymptotic expansion, while the latter will have the next order of smallness. We can now 
write 

qp{x) ~ res x^'^Ai{s), x ^ 0. 
Formula (31) allows the relation 

2 



'mm 



to be used in the calculations. 
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Turning to Eq. (29), we will obtain the final expression for the low-energy asymptotics 
of the function qp{x): 

qJx) ~ — — — , x^O. 34) 

Recall, however, that the function qp{x) is defined at x ^ 7inin ^ 1- Therefore, it is 
necessary to understand why Eq. (34) can be used in this case. The expression x ^ 
implies that the values of the argument under consideration arc small compared to 7^^. 
Recall that |-P(s)| < 1 in the regularity band. To qualitatively estimate ^tr we will set 
P{s) = in Eq. (29), which formally corresponds to ^ = 0. It follows from Eqs. (27) 
and (28) that the inverse transform of the derived function is proportional to the function 
(fi{x/-ftr), where 

koA 

Since the transition region for the function (p{y) is y ~ 1, the region of low energies is 
specified by the inequalities 7inin < 7|| < 7tr-- For a proper consideration, this region should 
not be vanishingly small. Therefore, it should be assumed that 

7min < Itr- (35) 

For typical values of ko ~ 10^, a ~ 100, and A ~ 10, we have a characteristic transition 
Lorentz factor •jtr ~ 10^. We see that criterion (35) definitely holds for a typical minimum 
Lorentz factor 7min ~ 100. 

As regards the high-energy asymptotics, we will not need it in the subsequent calculations. 
If necessary, it can be obtained by a method similar to that described above, but the path 
of integration should be transferred through the right bounding pole. Note only that the 
high-energy asymptotics is proportional to x~^'^'^'', where Smax is defined by Eq. (32). It has 
a power-law form rather than an exponential one, consistent with the conditions described 
after Eq. (24). 

Finally, let us calculate the required function (5(7min) (15). No direct integration of 
the function qp{x) is possible. However, the Mellin transform of the function Q{x) can be 
calculated using the relation 



M'"*(s)=9Jl[g(x)](s) = 



s 
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For the function Q{x), we then have 



Q{x) ^m-^[M''''{s)]{x). 

The function A1*"*(s) is meromorphic in the entire complex plane and regular in the band 
•^mhi < cr < SmL) where = Smin — 1 and = Smax — 1 are the left and right bounding 
poles, respectively. It remains to transfer the path of integration beyond the left bounding 
pole ^ 1 and to find that 

Q{x) ~ xqp{x), X — > 0. 

It makes sense to consider this asymptotic expression at x <^ ^ytr- Using criterion (35) and 
taking 7inin as x, we have 

Q(7min) = ^-4rk- (36) 

3 a m A 

It is easy to see that the source of electron-positron pairs Q depends linearly on the total 
intensity of the curvature radiation (16): 

Q - XW^urvB, (37) 

where 

X 3 A 1 

^= or^ — ^ 1- 

8 In A 7i„in 

In order of magnitude, A ~ 0.1. 

Above, wc considered the production of electron-positron pairs by curvature and syn- 
chrotron photons. The presence of synchrotron photons is important when considering an 
electron-positron cascade in the polar cap for a closed magnetosphere completely filled with 
plasma. In this case, even despite the electric field screening, apart from the primary parti- 
cles accelerated in the polar gap, there are two generations of secondary particles: the first 
generation of secondary particles produced by the curvature photons from primary particles 
and the second generation of particles produced by the synchrotron photons emitted through 
the transition of first-generation particles to the zeroth Landau level. In our case, we have a 
significant longitudinal electric field that continuously accelerates the newly produced par- 
ticles. For this reason, it makes no sense to consider the plasma production cascade as a set 
of successive generations, because the produced particles, whatever generation they formally 
belong to, immediately begin to play the role of primary particles after their acceleration. 
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Consequently, the presence of synchrotron photons in the case we consider is all the more 
important, because we formally have the case corresponding to an infinite number of par- 
ticle generations. In other words, the second term on the right-hand side of the integral 
equation (13) is much larger than the first one. From qualitative considerations, we can 
write the inequality 

Q » (38) 

with the characteristic number of curvature photons produced per unit time appearing on 
its right-hand side. Inequality (38) means that the existing synchrotron radiation makes a 
major contribution to the total pair production rate. Using Eqs. (37) and (38), we obtain 

7min « (39) 

In particular, given the inequality InA > 1, condition (35) follows from inequality (39). 
Since A ~ 10, inequalities (35) and (39) may be considered to be essentially equivalent. 
They both always hold, because it is implied that inequality (38) holds. Otherwise, in the 
presence of only curvature radiation, we should have formally set ^ = in Eq. (13), which 
would immediately give the solution. 

When the condition of quasi-stationary motion W^urv = E\\ is met, we have 

Q = A|E-B|. (40) 



3. THE INITIAL STAGE OF CASCADE DEVELOPMENT 

Let there be a primary electron-positron pair produced by a photon from the external 
cosmic gamma-ray background in a neutron star magnetosphere. In the presence of an 
external longitudinal electric field, the particles moving virtually along a magnetic field line 
begin to fiy apart in opposite directions and to undergo acceleration, reaching the Lorentz 
factor 7o. The emerging curvature radiation will be absorbed in the magnetic field through 
the production of secondary electron-positron pairs. We cannot use Eq. (36) from the initial 
time t — 0, because this requires that the particle number density changes little in the 
photon mean free path. Consequently, to be able to consider the electron-positron plasma 
generation processes using the hydrodynamic equations at least requires that a time equal 
in dimensionless units to the photon mean free path elapses since the initial time t — 0. 
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The particles of the electron-positron pair produced by this photon will then give secondary 
curvature and synchrotron radiation, which also produces pairs and the source of particles 
from the integral equation (13) can subsequently be used. 

Let us find the characteristic distance li from the production point of a primary electron- 
positron pair to the production point of the first secondary electron-positron pair. It can be 
represented as the sum of three parts, 

ll la -\- Ig -\- Icurvi 

where la = Ji/E is the distance at which the particle acquires some Lorentz factor 71 ^ 7||, 
le — p/71 is the formation length of curvature radiation, and Icurv = lf{kcurv) is the mean 
free path of the curvature photon emitted by the particle. Here, we introduced the mean 
free path 

^^^^^ ^ Ta 

dependent on the photon energy k [26]. The distance li is a function of the Lorentz fac- 
tor 7i. The minimum value should be chosen from the various /i, because efficient plasma 
multiplication begins at times t > minZi(7i). The extremality condition dli/d'yi — gives 

71 

the Lorentz factor 

'aE» 



and the distance 



1/4 

71 = (2p)V2 ( ^ ) (41) 



h = t^. (42) 



3 El 

Formulas (41) and (42) are valid under the following condition: 



^ < 1. 43 

a 

Since the expression on the left-hand side of inequality (43) at ~ 10~^, B ~ 10~^, and 
A ~ 10 is equal in order of magnitude to lO^'^ — 10^'^, the inequality itself always holds. Let 
us take p ~ 10^^. The Lorentz factor is then ji ^ 3 x 10^ and the corresponding distance is 
Zi 3 X 10^^, which corresponds to /i ~ 10 m in dimensional units. 

Next, we can consider the production of an electron-positron plasma in the hydrodynamic 
approximation by specifying linear electron and positron densities in the form 

Nr'-^^dih-\z\), (44) 
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where 9{x) is the theta function, as the initial conditions. It is important to note that the 
electron and positron densities as functions of the longitudinal coordinate z are determined 
to within distances of the order of the photon mean free path relative to the electron-positron 
pair production. In the case under consideration, the distance /i acts as this length and it 
does not matter at all which specific function N^^^^^z) is taken on such small scales. The only 
requirement is that the integration of this function over the longitudinal coordinate gives 
the total number of electron-positron pairs. In Eq. (44), we took a uniform distribution of 
particles in longitudinal coordinate. 

It was shown in [33] that in the interval from the production time of the first secondary 
pair Ti ~ li to the time of full particle acceleration r^t, only a few tens of electron-positron 
pairs are produced at i^y ~ 10~^ and a/A ~ 10. The corresponding ratio of the times is 
Ti/Tst ~ 0.9. After the time Tst, the particle ceases to accelerate and emits 

_ _5_ajo 

of curvature photons per unit time. This expression is derived by integrating Eq. (4) over 
all energies k from to oo. The characteristic mean free path is given by the expression 

_ 8 a 

''curv 9*^^** A' 

For a/A ~ 10, the ratio of the curvature photon mean free path to the particle accelera- 
tion time is Icurv/Tst ~ 0.06 <^ 1. Consequently, once the stationary Lorentz factor 70 has 
been reached, the production of electron-positron pairs with the source Qcurv begins almost 
immediately. In this case, the ratio of the mean free path to the formation length of the cur- 
vature radiation is large, ~ 10^; therefore, the radiation formation length may be neglected. 
Qcurv can be taken as the source of electron-positron pairs only at times Tgt < t < 2ti, 
because the contribution from the curvature photons produced by secondary particles will 
subsequently appear. 

Thus, the cascade process of electron-positron plasma multiplication will subsequently 
begin. In this case, the effective switch-on of curvature radiation can be assumed to occur 
Tst after the particle production if 7|| <^ 70- It is important to note that the minimum 
mean free path for synchrotron photons at A ~ 10 coincides (in dimensionless units) in 
order of magnitude with the time of full particle acceleration [33] . Hence, the production of 
secondary electron-positron pairs by the curvature and synchrotron photons emitted by the 
particle begins almost simultaneously. 
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4. THE TIME DELAY AND THE EFFECTIVE LOCAL SOURCE 



It can be concluded from the preceding section that the particles of each newly produced 
electron-positron pair contribute to the total pair production rate some delay time r after 
its production. We can then write the differential-difference equation 

^ = 2QiV.(t-x) (46) 

with the initial conditions 

where A^o is the number of electron-positron pairs at the initial time. Here, we denote the 
total number of electron-positron pairs by 



oo 

Nj:{t)= J N±{z,t)dz, 



where N+{z,t) and N^{z,t) are the linear positron and electron densities, respectively. 
Let Nt- — 2Qt be the number of particles produced in the time r by the photons emitted 
by one particle. It is also convenient to pass to the time t„ = t/r normalized to the delay 
time. Equation (46) can be solved using the Laplace transform (see, e.g., [34]): 

iVE(t)=iVo$^^(tn-p)^ t^T. (47) 

This exact solution allows the total number of electron-positron pairs to be found but does 
not allow the distribution of electrons and positrons in longitudinal coordinate to be found. 
In addition, it is applicable until the screening of the external longitudinal electric field 
begins. The hydrodynamic equations suggest the existence of a dependence of the pair 
production rate at some point z and at some instant of time t on the electron and positron 
densities taken at the same point and at the same instant of time. In other words, the 
equations under consideration arc local in space and time. Nevertheless, there exist both 
spatial and temporal separations of the photon production and absorption points. If the 
change in electron and positron densities on characteristic photon mean free paths is small, 
then the spatial nonlocality condition imposes no significant constraints on the applicability 
of the hydrodynamic equations, except the impossibility of considering the plasma generation 
on short time scales (see the preceding section). In contrast, the temporal nonlocality 
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condition is significant, because, as we will see below, a large number of photons are produced 
in the delay time: 

> 1. 

Consequently, disregarding the temporal nonlocality would cause the electron-positron pair 
production rate to be grossly overestimated. 

We need to find some effective local source Q^-^^ that can be used in the hydrodynamic 
equations. To calculate Q^-^^, we will assume that all of the parameters characterizing the 
plasma tube change little at a distance of the order of the photon mean free path If. We 
can then select some tube segment with a length equal to several // in which the local 
homogeneity condition is met. For the tube segment under consideration, the source Q 
and the linear electron and positron densities N± are essentially invariable as functions of 
the longitudinal coordinate z. For the total number of electron-positron pairs in the tube 
segment under consideration, we can then write Eq. (46). 

Let us investigate solution (47). Let us introduce a function 

5(p)=pln(e7V,e(p)), 

where 

«!>) = 7-1- 

Replacing the summation by the integration over the variable p in Eq. (47) leads to an 
asymptotic estimate of the series 

N^^it) ^NoJ exp(5(p)) dp. (48) 



Let po be the point of maximum of the function S{p), i.e., 

SiPo) = max S{p). 

The equation 

i_expi = 7V„ (49) 

where — C(Po)) then follows from the extremality condition S'{po) = 0. At A^^ ^ 1, the 
inequality <^ 1 holds and can be easily found by the method of successive approxima- 
tions: 

^0 ^ [lnA"^-ln(lnA"^)]-^ 
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We then have the point of maximum 

t. 



n 

Po 



1 + ^0 

We can now estimate the integral in Eq. (48) by the Laplace method [35] to obtain 

N^'{t) = ^{eN^Cor. (50) 

The figure shows the time dependence of {N^^{t) - A^s(i))/-/Vs(i) for = 10^ We see 
that our asymptotic estimate under the condition i„ ^ 1 closely approaches the exact 
solution, with the error tending to zero as t increases. At lower A^^-, the relative error 
only decreases. The function iVs(t) is only continuous rather than infinitely differentiable. 
This is due to the stepwise "switch-on" of the next generation of particles contributing to 
the total electron-positron pair production rate at time intervals r. The most important 
thing here is that the plot of N^^{t) intersects the plot of A^s(^) in an arbitrarily chosen 
interval nr < t < {n + 1)t, n — 0, 1, 2, ... . Consequently, the derived smooth monotonic 
approximation N^{t) smoothes NY,{t). Even differing from the exact solution in certain 
time intervals, it properly describes the mean rate of increase in the number of particles and 
allows the error accumulation in the number of produced pairs on long time scales to be 
avoided. 

Differentiating Eq. (50) with respect to time using relation (49) gives 
where we specified the effective source of electron-positron pairs 

Thus, the number of electron-positron pairs will increase exponentially with the effective 
local source Q^^^ in each selected segment of the plasma tube. 

It remains to determine the delay time r and the initial source Q. Consider one of the 
particles of the produced electron-positron pair. If we exclude the synchrotron photons from 
consideration, then the curvature photons emitted by the particle will begin to contribute 
to the production of secondary electron-positron pairs after the time Tgt of full particle 
acceleration. In this case, Qcurv defined by Eq. (45) should be taken as the source Q. This 
case corresponds to the presence of a strong magnetic field B > 1. If, alternatively, we 
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also include the synchrotron photons in consideration, then, as follows from the preceding 
section, they will begin to contribute to the pair production on time scales of the same order 
of magnitude. In addition, when deriving the integral equation (13), we assumed that each 
newly produced particle accelerated to the Lorentz factor 70 and began to act as the source of 
primary curvature photons in the expression for qcm-v This will also occur the time Tgt after 
the particle production. Before this instant of time, the particle under consideration is not 
involved in the production of secondary electron-positron pairs and, consequently, makes no 
contribution in the form of the fraction of the total pair production rate per particle defined 
by (36) and Eq. (40) following from it. Thus, Eq. (40) and the time Tgt should be taken as 
the source Q and the delay time r, respectively. This case corresponds to the presence of a 
weak magnetic field 5 -C 1. We then have 

Nr., = 2A70S > 1. 

At S ~ 0.01 and 70 ~ 10^, we obtain N^^^ ~ 10^, i.e., the inequality N^^^ > 1 holds. The 
corresponding value of ^0 ~ 0.1. Note that when the electric field is subsequently screened, 
this inequality will also hold even if the longitudinal Lorentz factor gradually decreases down 
to the small value 70 ~ 10^. The required effective local source will be defined by Eq. (51), 
in which r = Tgt and ^0 is the solution of Eq. (49) for = N^.^^ . Finally, we can write 

Q'^' - (52) 

5. CONCLUSIONS 

The absorption of a high-energy photon from the external cosmic gamma-ray background 
in the inner neutron star magnetosphere triggers the generation of a secondary electron- 
positron plasma. The existence of a strong longitudinal electric field is a necessary condition 
for the development of a cascade. Allowance for the pair production by curvature and 
synchrotron photons in a weak (compared to the Schwinger one) magnetic field necessitates 
considering the integral equation (13). The contribution from synchrotron photons to the 
source of electron-positron pairs in the case of a weak magnetic field is decisive. In this case, 
the contribution from bremsstrahlung photons may be neglected. The integral equation is 
solved using the Mellin transform. The transfer of the path of integration beyond the pole 
bounding the analyticity band of the transform of the solution on the left and the method 
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of residues allow one to obtain the low-energy asymptotics of the solution and to show its 
applicability for typical conditions in the neutron star magnetosphere. If the conditions of 
quasi-stationary motion are met, then the solution is reduced to Eq. (40). In the case of a 
strong magnetic field B > 1, because of the absence of synchrotron photons, the source of 
electron-positron pairs is determined by the intensity of the curvature radiation and has the 
form (45). Immediately after the absorption of the primary gamma-ray photon triggering a 
cascade in the magnetosphere, the hydrodynamic equations are inapplicable. The particles 
of the produced primary electron-positron pair begin to move in opposite directions almost 
along a magnetic field line and reach a Lorentz factor of ~ 10^ in the initially vacuum electric 
field. The charged particle acceleration time, the curvature radiation formation length, and 
the photon mean free path relative to the electron-positron pair production are finite. This 
determines the production of secondary pairs at a distance of ~ 10 m from the primary 
pair. This distance is comparable to the distance of particle acceleration to a stationary 
Lorentz factor. Subsequently, the next generation stage begins, when the pair production 
occurs due to intense absorption of the curvature photons produced by the particles of the 
primary pair during their acceleration. Since the production rate of curvature photons by 
a fully accelerated particle is high, the effective switch-on of curvature radiation occurs 
Tst after the electron-positron pair production. The same is also true of the synchrotron 
photons: their minimum mean free path is comparable to the distance at which the particle 
reaches a stationary Lorentz factor. The absorption of the synchrotron photons produced by 
the particles of secondary electron-positron pairs when passing to the zeroth Landau level 
gives rise to the main stage of electron-positron plasma generation when its rate is highest. 
The high secondary plasma generation rate due to the existence of a strong accelerating 
electric field necessitates allowance for the time delay r. After this time, the particles of 
each newly produced electron-positron pair begin to contribute to the total pair production 
rate through the absorption of the curvature and synchrotron photons emitted by them in 
the magnetosphere. Solving the differential-difference equation (46) allows one to calculate 
the effective local source of electron-positron pairs (52), using which the time delay can 
be properly taken into account. Significantly, the delay leads to a great reduction in the 
electron-positron plasma generation rate compared to the case of zero time delay. The 
effective local source (52) is much weaker than both the initial source (45) for S > 1, when 
there is no synchrotron radiation, and the initial source (40) for S <^ 1, when the intensity 
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of the synchrotron radiation exceeds that of the curvature one. In this case, the delay time 
T for both weak and strong magnetic fields is almost the same and is comparable in order 
of magnitude to the particle transition time Tgt to a quasi-stationary regime of motion. 

A detailed study of the formation of a plasma tube in a neutron star magnetosphere 
through the cascade multiplication of an electron-positron plasma triggered by the absorp- 
tion of a high-energy photon from the external cosmic gamma-ray background with allowance 
made for the screening of the external electric field requires a separate consideration. 
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